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Abstract. We propose that Kozai's phenomenon is responsible for the long-term evolution of stellar orbits near a supermassive 
black hole. We pursue the idea that this process may be driven by a fossil accretion disc in the centre of our Galaxy, causing 
the gradual orbital decay of stellar trajectories, while setting some stars on highly elliptic orbits. We evolve model orbits 
that undergo repetitive transitions across the disc over the period of » 10' years. We assume that the disc mass is small 
compared to the central black hole, and its gravitational field comparatively weak, yet non-zero, and we set the present values 
of orbital parameters of the model star consistent with those reported for the S2 star in Sagittarius A*. We show how a model 
trajectory decays and circularizes, but at some point the mean eccentricity is substantially increased by Kozai's resonance. In 
consequence the orbital decay of highly eccentric orbits is accelerated. A combination of an axially symmetric gravitational 
field and dissipative environment can provide a mechanism explaining the origin of stars on highly eccentric orbits tightly bound 
to the central black hole. In the context of other S-stars, we can conclude that an acceptable mass of the disc (i.e., < 1 per 
cent of the black hole mass) is compatible with their surprisingly young age and small pericentre distances, provided these stars 
were formed at r < 10^ gravitational radii. 
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1. Introduction 

' A dense nuclear cluster surrounds the centre of our Galaxy. 
Tracking the rapid proper motion of individual stars in the cen- 
tral arcsecond (i.e., about 0.04 pc off Sgr A* object) provides an 

. essential tool to determine the central mass M, 3.5 x IO^Mq 
(Genzel et al.'2003', Ghez et al. '2003a). So far, the motion of S2 

. (SO-2) has been measured with the best accuracy over a suffi- 

' cient duration of time (Schodel et al. l2002ll2003l . It is a highly 
eccentric, e - 0.87 Keplerian orbit with the orbital period of 
15.2 years and an apocentre distance 0.01 pc. Other newly dis- 
covered stars strengthen the case for a supermassive black hole 
in Sgr A*, and for the surprisingly young stellar population 
around it. The rather high mass and low age present a chal- 
lenge to star-formation theories (Milosavljevic & Loeb 2004 
Eckart et al.'2004'). In this paper we discuss one of the possible 
mechanisms that could help to resolve the issue. 

S2 is a member of the central Galactic cluster for which ra- 
dial velocity as well as proper motion have been reported by 
Ghez et al. 1 2003b l. Combination of the two pieces of informa- 
tion allowed these authors to determine the orbital parameters 
from the full three-dimensional motion and to break the am- 
biguity in the inclination angle. Trajectories and spectroscopic 
measurements of other stars in this region have not yet reached 
comparable accuracy, however, there are good prospects for 
future progress and strong indications of the puzzling situa- 



tion: a number of S-stars are on close orbits near the black 
hole but they seem to be of young age, contrary to expec- 
tations. Several mechanisms have been proposed that could 
bring stars to the neighbourhood of the central black hole (e.g., 
Hansen & Milosavljevic '2003 ; McMillan & Portegies Zwart 
2003, Alexander & LivioJ004 1. In spite of this effort the origin 
of close stars is not well understood. Models have difficulties in 
reconciling different aspects of the Galaxy centre - on the one 
hand it is a low level of present activity, indicating a very small 
accretion rate, and on the other hand the spectral classification 
suggests that these stars have been formed relatively recently. 
While in-situ star formation is problematic because of strong 
tidal forces operating near the central black hole, two-body 
gravitational relaxation could in principle set stars on plung- 
ing orbits. Nonetheless, the latter mechanism does not seem 
to be efficient enough because a rough estimate of the relax- 
ation time-scale, f, > 1 Gyr, is orders of magnitude longer than 
the estimated age of the S-stars (x 10 Myr). We neglect the 
two-body gravitational scattering and the effect of gravitational 
relaxation in the present analysis (see Alexander & Hopman 
p003 for a more elaborate discussion of gravitational relaxation 
in the Galactic Centre). 

The interaction between stars and a fossil accretion disc 
is a cumulative effect which provides a viable mechanism of 
transporting stars towards the black hole, while stripping their 
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outer atmosphere away. These two effects are mutually linked 
and they complement one another: whereas the stellar orbits are 
evolved by the interaction with the disc, both gravitational and 
hydrodynamical as discussed in this paper, the stars themselves 
as well as the composition of the gaseous environment are also 
affected (Artymowicz et al. 1993, Armitage et al. 1996 1. The 
main appeal of this scheme lies in the fact that it naturally com- 
bines several influences which are needed in order to explain 
the efficient migration of stars to the centre: the proposed sce- 
nario sets some stars from the sample on elliptic orbits for a 
fraction of their lifetime (the gravitational effect of the disc), it 
also causes continuous migration (hydrodynamical effect of the 
orbital dissipation), and finally it affects stars by removing their 
atmospheres. All these effects are determined by the assumed 
density profile of the disc and they are thus interconnected. 
Furthermore, it was proposed that S-stars have undergone re- 
juvenation by the interaction with the environment or by tidal 
squeezing during close encounters with the black hole, thereby 
hampering the spectral classification and mimicking young age 
(e.g. Morris[R93 Alexander & Livio 120(311 . Notice, however, 
that we neglect these peculiarities, nor do we take into account 
the likely effect of the disc itself evolving. 

Various authors have investigated the long-term orbital evo- 
lution of a stellar satellite that is either periodically crash- 
ing on the black-hole accretion disc (e.g., Syer et al. .1991. 
Vokrouhlicky & Karas IT9931 I1998al Rauch [T9951 Subr & 
Karas 1999 1 or stays embedded in the disc plane (Artymowicz 
[TW3. Hameury et al. .l994 Karas & Subr 2001 1. The above- 
mentioned papers address the star-disc interaction in active 
galactic nuclei, whose accretion discs are supposed to be rather 
dense, however, it has been suggested (Levin & Beloborodov 
120031 Nayakshin & Sunyaevl2()Q3J that also our Galaxy centre 
contained a similar gaseous disc at some stage of its history. 
Some form of an inactive disc or a dusty torus may still ex- 
ist as a remnant of past activity. If such a disc existed in the 
Galaxy centre, it could capture passing stars from unbound or- 
bits, or it could directly supply the medium where new stars 
form and evolve while sinking to the black hole. This mech- 
anism has been explored under conditions relevant for active 
galactic nuclei and quasars (Collin & Zahn .1999. Goodman 
& Tan l2004t and for Sgr A* (Genzel et al. l2003l . Likewise 
it has been speculated that star-disc collisions may contribute 
to the featureless X-ray variability by creating bright spots 
on the disc surface and fountains of ejected material above it 
(Zentsova'1983' Karas & Vokrouhlicky 1994, Vilkoviskij & 
Czernv .2002) . Cuadra et al. (2003 1 and Nayakshin et al. (2 003 
I2004bt have pointed out that a similar mechanism may be re- 
sponsible for X-ray flares, the occurrence of which was first 
reported by Baganoff et al. (2001 1. A sufficiently long observa- 
tion could reveal periodicities of the stellar motion connected 
to precession of the orbit, and from them one would be able to 
infer the angular momentum of the central black hole (Semerak 
et al. [T999l Aschenbach et al. l2004l . 

We apply the model of star-disc collisions to the orbital 
evolution of the star S2, assuming a dissipative perturbation 
of its motion by crossing an inactive gaseous slab. We em- 
ploy a similar formulation of the model as described by Subr 
& Karas (.1999 J and Karas & Subr (.2001 J. The relative mass 



of the disc, e = M^jM,, is a free parameter for us. We 
adopt a plausible and quite conservative value of e in the range 
IQ-^ <, e < IQ-^ (see also Nayakshin & Cuadra l2004al Subr 
et al.'2004). Simulations are initiaHzed at n > 3 x IQ'^Rg (i.e. 
above the present S2 distance; = GM.jc^ 1.9 x lO^^pc) 
and we follow the sinking satellite as it proceeds to the centre. 
The gas-assisted drag tends to circularize the orbit and align it 
with the disc plane, while the Kozai ( 1962) mechanism drives 
eccentricity to a large value, characteristic of the present stage 
of S2. We employ this mechanism to explain the origin of the 
large orbital eccentricity and small pericentre. The assumed 
mass of gas distributed in the disc and the resulting perturba- 
tion of the predicted orbits are small enough to be consistent 
with Keplerian motion over one (or several) revolutions, which 
is the time span accessible to direct observational checks (see 
Mouawad et al. 2004 for an upper estimate of the mass outside 
the black hole). Clearly visible departures from perfectly ellip- 
tical orbits are expected on a substantially longer time-scale, 
K 10^ orbits. 

2. The orbit evolution 

Let the gravitational field of the central black hole be described 
in the Newtonian approximation and the disc be planar and ge- 
ometrically thin. Passages of a star through the disc last only a 
small fraction of the orbital period at the corresponding radius, 
and so they can be considered as instantaneous events at which 
the passing body exchanges momentum with the material ly- 
ing along its trajectory. Repetitive events cause the orbit to de- 
cay. Outside the disc plane the star moves freely and its orbit is 
found numerically by integration in the combined gravitational 
field of the central mass M, and the disc mass M^. This requires 
one to specify the surface density distribution as a function of 
radius in the disc and we assumed two cases for simplicity. 
First we employ a homogeneous slab Ed(^) - const. Its grav- 
itational field can be then expressed analytically (e.g.. Lass & 
Blitzer ,1983. Pierens & Hure .2004.) . which helps to speed up 
computations. Another variant, l^diR) oc -\/R(Rout - R), was 
used to mimic the surface density and the gravitational field 
of an annulus (Rin<R<Rom)- Given the value of Md, the non- 
sphericity of the gravitational potential will be smaller and its 
impact on stellar orbits weaker if the disc matter is concen- 
trated more towards the centre. However, we checked that a 
very similar effect to that described here is obtained even with 
the standard Shakura-Sunyaev density profile, i.e. decreasing 
with radius as I^diR) R^^^^. Most likely the actual density 
distribution will be somewhat different because the disc will be 
modified by self-gravity, stellar transitions and other effects.' 

The assumptions about the form of the disc are necessary 
for definiteness of examples, but they do not affect the nature 
of our results. If the trajectory is oriented in such a way that it 
crosses the disc repetitively, we observe the gradual dissipation 

' By the effects of self-gravity we understand the gravitational in- 
fluence of the disc on the stars. As mentioned above, we do not modify 
the internal structure of the disc itself although this may also be inter- 
esting to examine because self-gravity is expected to cause a clumpy 
structure of the disc, thereby reducing the efficiency of star-disc colli- 
sions. 
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of orbital energy and of angular momentum. As time proceeds, 
the orbit becomes circular and declined in the disc plane, but 
this overall trend is overlayed with oscillations of eccentricity 
and inclination. Occasional jumps toward high eccentricity oc- 
cur due to self-gravity of the disc, acting as a perturbation in 
our model. 



Kozai's method was developed to describe the evolution of 
mean orbital elements of a small body disturbed by the pres- 
ence of another mass in the system. Our aim here is to consider 
a restricted problem in which the origin of perturbation is an 
accretion disc and the star represents a test body (a hierarchical 
triple system; see Lidov & ZigUn ll975l Thomas & Morbidelli 
[1996 1 Kinoshita & NagaifT999'). Kozai's resonance is a partic- 
ular type of secular resonance occurring when the pericentre 
argument precession vanishes, w = for the test body. In or- 
der to disregard various short-time effects acting on dynamical 
time-scale, the method invokes averaging technique which is 
based on canonical transformation (Arnold ^989; Brower & 



Clemence 119611 Morbidelh I2002t . Although this formalism 
was suited and intensely applied to the problem of motion in 
the Solar system, it has met numerous applications elsewhere. 
For example, Holman et al. ( 1997 1 proposed that large eccen- 
tricities of the orbits of some extrasolar planets could be under- 
stood in terms of Kozai oscillations, which are caused by the 
disturbing influence of the secondary star. A similar approach 
was employed in the context of stellar dynamics around mas- 
sive black holes in galactic nuclei (Sridhar & Touma 1999 1 and 
in globular clusters (Wen .20033 . Blaes et al. (.2002) proposed 
that oscillations of eccentricity could be important for deter- 
mining the rate of binary black hole merging and gravitational 
wave emission expected from galactic centres. 

An important point to notice is the presence of dissipa- 
tion due to stars crossing the disc. The star experiences or- 
bital decay due to repetitive transitions changing the trajectory 
and pushing it across the separatrix. We consider two effects 
which both depend on the density of the gaseous environment: 
the energy dissipation is roughly proportional to the density 
of the medium and likewise the gravitational perturbation in- 
creases its impact proportionally to the mass contained in the 
disc. These two processes are therefore mutually connected. 
As far as star-disc collisions are concerned, relevance of the 
Kozai-type mechanism coupled with dissipation was originally 
discussed by Vokrouhlicky & Karas ( 1998b) and we refer the 
reader to their paper for further details. Here we briefly sum- 
marize the adopted approach. 



The origin of the disturbing gravitational influence is self- 
gravity of the accretion disc, hence the perturbation potential 
is already assumed axiaUy symmetric. Averaging is performed 
over the star ellipse. We denote the dominant potential of the 
central source as Vdr) = -GM,/r (the central black hole) and 
the disturbing component Vd{R,z) (the disc). Test-particle mo- 
tion in the referential background field is described by an un- 



terms of semi-axis a. Eccentricity and inclination are related 

by2 



T] COS / = c, with = Vl - 



(1) 



and c secularly constant (Kozai's integral). Along c - const 
orbital oscillations take place on substantially longer intervals 
than the dynamical period at corresponding distances, defin- 
ing the medium time-scale fs of the problem. A fairly accurate 
estimation of can be obtained by replacing our disc with a 
narrow ring of radius Rd and the same mass (e.g. Kiseleva 
etal. ■19981 . 
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where fp(fl) is the Keplerian orbital period and /^d— >^out char- 
acterizes the size of the disc. More precisely, c as well as a 
would have been conserved if the energy dissipation had been 
negligible. However, the orbital energy is gradually dissipated 
in our system and this ensures that the star reaches the separa- 
trix and passes through it, i.e., the orbital elements change their 
values and they eventually acquire a suitable combination. The 
overall orbital decay introduces the third time-scale, fc, and we 
assume that it is the longest time-scale relevant for our problem 
(this assumption imposes an upper limit on density of the disc 
but one can check that the constraint is safely met under astro- 
physically realistic conditions). The total duration is obtained 
by integrating the inverse rate of orbital decay from the outer 
edge of the disc down to (say) tidal radius. 



tc oc 
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where a < depends on the surface density profile. Its magni- 
tude is roughly oc 2(7?), but details are model dependent (Syer 
et al. 1991] Zurek et al. 199_fJ. Following the arguments in 
Karas & Subr J2001I) we estimate 
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Here we introduced £3 = E,/(0.32q) with 2» = Mt/S , and 
St= ttRI (column density, mass, cross-sectional area and ra- 
dius of the star, and analogically for the sun). S3 1 for S2, 
i.e. a fiducial value for young, massive stars. Clearly the pe- 
riod f(, is inversely proportional to rI and found of the order of 
fc ~ 10" Myr for typical values of the model parameters. 

As mentioned above, the main difference between the di- 
rect numerical integration of stellar orbits and the analyti- 
cal approach of the averaging technique lies in the way how 
the rapidly changing (fast) variables are treated. The mean 
anomaly plays the role of a fast variable in our problem. 



perturbed Hamiltonian, which is simply 'He 



iGM.fl-i in 



- In our notation, ; = 0° corresponds to an orbit co-rotating with the 
disc, while ; = 180° is for counter-rotating one. 
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Fig. 1. Contour lines of the disturbing function Vi(p, a>) (shown with equidistant steps between contour levels). The polar plot 
is shown with zero eccentricity at the origin p - Q (e = 0) and polar angle oj indicated on margins. In the case of negligible 
dissipation the orbital evolution proceeds along contours. In different regions of the graph oj either circulates over the full circle 
or it librates within a limited interval. A separatrix (thick line) forms a border between these regions. Four topologically different 
cases can be distinguished by the structure of contours, shape of the separatrix and the coiTesponding value of Kozai constant 
c: (a) c = (e < 1); (b) c = 0.2 (e < 0.98); (c) c = 0.8 (e < 0.60); (d) c = 0.9 (e < 0.44). These define the range of orbital 
oscillations of eccentricity and inclination according to eq. Q. Orbital energy has been set identical and equal to -GM,/{2a) 
with semi-axis a - 2.4 x lO*Rg in all panels. While on the medium time-scale the orbits remain attached to underlying contours, 
allowing us to estimate the range of variations of eccentricity, on the long time-scale the orbits are driven across contours. This 
gradual diffusion leads to topology changes at some moments. 



therefore, the aim is to find a canonical transformation under 
which the transformed Hamiltonian becomes independent of 
this quantity. The transformation depends on the form of the 
disturbing potential VdiR,z), and so it cannot be given explic- 
itly for an arbitrary profile of density in the disc. One can how- 
ever introduce a formal development with respect to the small 
parameter e <& I. After transformation and expansion to the 
first order we find 

Va^Vi^Vi+0(e). (5) 
The first term of the series can be written explicitly, 

= r (-fvdR,z)dv, (6) 
Ittt] Jo ^a/ 



where integration is performed over the true anomaly v of the 
star for one revolution around the black hole. At this level of 
approximation the mean semimajor axis stays constant and dif- 
fers from the osculating semimajor axis by short-period terms 
only (we therefore do not distinguish them in our notation). The 
problem is reduced to the evolution of the mean eccentricity e 
and pericentre argument a>, which are further constrained by 
Vi (e, a>; a, c) - const. When evaluating this condition, a and c 
are parameters. 

Levels of Vi - const are shown in Figure [2 We use non- 
singular Poincare coordinates k(e, a>) and h(e, w), defined by 
+ = 2(1-7]) and h/k - tanw. In order to draw these 
contours we construct a polar graph in (p, co) with the azimuthal 
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Fig. 2. Time evolution of a representative trajectory, which evolves according to the mutual interplay of gravitational and dis- 
sipative interaction with an accretion disc. Mean orbital elements are shown in different panels. The semi-axis a(T) (solid line) 
monotonously decays as the star sinks to the centre. Inclination /(t) and eccentricity e(T) (plotted with points) oscillate with rather 
large amplitudes. We also show a reference orbit (dashed line), which was computed by neglecting the gravitational field of the 
disc. The oscillations and the Kozai resonance are due to the disc gravity; hence they affect only the exact trajectory and not 
the referential one. The orbital decay is due to dissipation, hence it affects both solutions. Both cases start with the same initial 
values: a; - 3.5 x IQ*R„, - 0.35 and /j = 84°. Zero time has been set at the moment when the orbital parameters of the exact 
trajectory change their course abruptly, attaining the values consistent with those reported for the S2 star. We also show the graph 
of inclination vs. eccentricity (bottom right) with an elliptical arc segment corresponding to Kozai's c - 0.275. The transition 
to a different regime of Kozai's oscillations occurs along this segment and a detailed subplot of the indicated rectangular area is 
enlarged in Fig.|3l 



angle < oj < 2n, polar radius p = y/h^ + and the origin 
p - coinciding with the circular orbit (the mean eccentricity 
increases with polar radius, e cc p near origin). The disc surface 
density was assumed constant. 

Rigorous proof of convergence of the expansion (|5j and 
an estimation of accuracy of the first term Vi in the series is 
a difficult task which cannot be assessed for a general pertur- 
bation potential. However, a partial check can be achieved by 
comparing the analytical result with numerical integration of 
exemplary trajectories (we will show such a comparison in the 
next section). Typically, results based on a truncated part of the 



series are valid over a limited time span and they fail to predict 
the system evolution for infinite duration. 

We conclude this section by summarizing the main points 
in our perspective of how Kozai's phenomenon could help to 
solve the problem of S-stars in Sgr A*: 



(i) Axially symmetric (non-spherical) perturbation of the cen- 
tral gravitational field is the essential element of the pro- 
posed scheme. We attribute its origin to self-gravity of the 
gaseous environment with the geometry of the disc. 
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Fig. 3. Detail of the trajectory from Fig. |2] showing that the 
orbit (solid line) oscillates almost precisely along the lines of 
Kozai's integral (several contours are plotted with dashed lines 
and the corresponding values of c = const are given). This in- 
dicates that the disc mass is small but non-zero, and its gravity 
provides a perturbation of the central field of the black hole. 
The orbit diffuses across the c-curves due to orbital decay by 
dissipation. 

(ii) Stars of the nuclear cluster follow nearly Keplerian orbits 
around the central black hole. Free motion in the central 
field is disturbed by the disc gravity. 

(iii) Star-disc collisions lead to dissipation and gradual 
sinking toward the black hole, the decay rate being 
faster for counter-rotating and highly eccentric trajectories. 
Eccentricity of highly inclined orbits is pumped to large 
values via Kozai's mechanism. Consequently, the pericen- 
tre distance shrinks. 

3. Results 

Examples of integration of the evolving trajectory are shown 
in Figure |2] where the orbital oscillations and the progres- 
sive sinking of the star to the centre are evident. We introduce 
rescaled time t as a convenient measure, related to time t in 
physical units by 



M. 



3.5 X 10''M( 



O 



(7) 



The characteristic time-scale of the oscillations depends on 
Md/M, ratio and is generally much longer that the dynamical 
time at the corresponding radius (fs is of the order of thousands 
of orbital revolutions in our case). On the other hand, the rate 
of orbital decay is proportional to S» as indicated by eq. @, 
and so our results can be readily scaled to different stellar types. 
Furthermore, scaling can be performed to different values of 
M, provided that the ratio Mi/M, is retained, as well as the 
disc dimensions and the orbital semimajor axis in units of Rg. 

Now we set the inner edge of the disc at = 6R„ and the 
outer edge 7?out = 5 X 10'^ R„. The disc mass is 
which for M. = 3.5 x IO^Mq gives Ed = 2.9 x lO^gcm^^ ^j. 
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Fig. 4. The probability p{t; e^) of attaining eccentricity higher 
than ec for the same example orbit as in Figure|2l Three curves 
are shown for different values of ec- The eccentricity reaches 
large values during the period starting at t = and lasting 
about 6 Myrs, and it stays above = 0.7 for a considerable 
fraction of time. 



surface density. In this case, Fig.|2demonstrates how eccentric- 
ity is pumped to large values and, consequently, the orbital de- 
cay is accelerated. The semimajor axis decreases by a factor of 
5 in several million years, even for a relatively moderate choice 
of the disc mass and density. The effect is thus potentially im- 
portant, as it is able to bring stars near the black hole and put 
them on eccentric orbits. In these computations we employed a 
fourth-order Runge-Kutta integrator with adaptive stepsize. In 
order to reduce the duration of runs we pre-generated the field 
components in a Cartesian grid covering a rectangular region 
of 1 pc X 1 pc. In some regions, close to the disc, the field lines 
are highly distorted. Therefore, we used additional sub-grids 
that provide fine resolution where needed. Field values in be- 
tween mesh points were obtained by bilinear interpolation. As 
detailed below, we can be fairly confident of sufficient accuracy 
of the numerical computation because the resulting trajectories 
agree with expectations based on analytical arguments of the 
perturbation method. 

The orbital oscillations, which we observe in graphs, call 
for some kind of probabilistic interpretation of the orbital pa- 
rameters. This can be inferred from non-uniform distribution 
of dots which represent fluctuations of e(T) and ;(t) in Fig.|2] 
Each dot represents the position in the phase space at a corre- 
sponding moment of time, and the density of dots determines 
the probability that the star occurs with given values of param- 
eters.^ From Figure|3and especially from the detail in Figure^ 
we see that the oscillations of eccentricity and inclination are 
bound together by Kozai's quasi-integral c = Vl-? cos / - 
const. 



To avoid fake structures in the graph, dots are plotted with random 
time steps that are distributed uniformly over (0, t^.^,^), where t^^^ is 
an interval longer than the orbital time and shorter than characteristic 
time f. of orbital oscillations. 
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Fig. 5. The contour plot of disturbing potential Vi(p, ai) - const shown at two moments of time. The exact (numerically inte- 
grated) trajectory is also plotted (by points). The left panel corresponds to the moment just before the numerical trajectory crosses 
the separatrix: r - -O.lSMyr, a = 2.50 x IQ'^Rg. The corresponding value of Kozai's c = 0.2733 and the maximum eccentricity 
Sm-dx - 0.962. On the other hand, the right panel captures the orbit after crossing the separatrix: t = O.lSMyr, a = 2.46 x lO'^Rg, 
c = 0.2795, emax = 0.960. 



Another representation of the orbital oscillations is given 
in Figure |4] where we plot the time evolution of probability 
p(T;ec) that eccentricity has exceeded a given threshold e^. 
This graph demonstrates how a sinking star can, indeed, be 
set on a highly elliptic trajectory at certain stages, while grad- 
ually approaching the black hole. Fig. |2 helps to answer the 
question of what fraction of time the orbit actually spends 
at the stage of high eccentricity. For example, by integrating 
Jp(t; - 0.8) dr over the first 6 Myrs we find the star spends 
approximately half its time with an eccentricity larger than 0.8. 
Then the orbit becomes progressively circularized and the star 
descends below IQ'^R^. 

b 

Contour levels of the disturbing potential are shown in 
Figure|5l We construct this graph using the same procedure as 
described for Fig.^ Two panels correspond to slightly differ- 
ent moments; notice that there is a small change of the contour 
line structure because of continuous dissipation. In these plots 
we also include the exact (numerically integrated) trajectory. 
Clearly the adopted approximation is sufficiently accurate: the 
numerical trajectory follows its underlying contour at each cor- 
responding panel, as it should. An adiabatic diffusion across 
Vi - const lines is slow and it cannot be resolved in this figure, 
where cj circulates over the whole interval (0, In) on the time- 
scale of the order of several hundred orbital periods. However, 
there is an evident difference between the two cases: in the left 
panel the trajectory oscillates near the origin keeping small ec- 
centricity, while on the right it has crossed the separatrix. As a 
result of this change, the coiTesponding orbital energy is also 
slightly different, whereas the mean eccentricity starts to vary 
over substantially larger intervals. This also explains how the 
averaged model can be used to estimate the maximum eccen- 
tricity attained by the star. 

In the above example the Kozai oscillations were only used 
to set the S2 orbit into a highly eccentric state. Indeed, it is 



evident from the orbit evolution at t > that its decay is con- 
siderably accelerated by this process. This is due to the fact 
that a star on a highly eccentric orbit passes through the disc 
with large relative velocity. The same argument suggests that 
counter-rotating stars suffer more from the orbital dissipation 
than co-rotating ones. Moreover, the efficiency of dissipation 
due to large inclination is strengthened because the separa- 
trix prevents the overturning of counter-rotating orbits to co- 
rotation. 

Hence, we found it convenient to explore the possibility 
that the S-stars have been formed on nearly circular orbits far- 
ther form the centre. Then, due to the simultaneous energy dis- 
sipation and the influence of the perturbing gravitational field 
the stars undergo orbital decay. Sample trajectories in Figure|6l 
were integrated assuming Rout -2x 10^ Rg and M^jM, - 0.01. 
Two different lines represent the orbital decay of a star interact- 
ing with different Ed(^) profiles starting with identical initial 
values of orbital elements. 

The region of efficient star formation is usually placed at 
r X 10^ Rg, however, the minimum radius is not known very 
precisely. Collin & Zahn H^Wi argue, using the model of a 
marginally unstable self-gravitating accretion disc, that the star 
forming process can take place at distances > 10^/?g. Although 
their work was intended for quasars, it indicates that the mech- 
anism discussed here could be relevant also for Sgr A* if there 
is enough material dispersed around the black hole, or if it was 
there in the past. Star forming process occurring at roughly 
IQ^Rg was also considered by Nayakshin & Cuadra (.2004a) . 
Furthermore, it was suggested (Genzel et aI. I2003t that young 
He I stars rotating in two coherent rings of radii > 10^ Rg could 
have formed during a recent star forming episode triggered by 
a collision of different gaseous streams. Our two examples in 
Figure|6lindicate that the exact form of the gas distribution play 
an essential role, both for the orbital decay and for oscillations. 
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Fig. 6. Left: The orbital decay due to interaction with the disc. Two cases are shown corresponding to a homogeneous disc 
2d cx const (solid line), and oc -\JR(Rout - R) (dashed line) representing an annulus. The rate of gas-assisted orbital decay 
speeds up at the moment of crossing the separatrix (time r = 25 Myrs and r = 47 Myrs, respectively). Right: the probability for 
these orbits to have eccentricity larger than 0.8. Both trajectories start with identical initial values, Ui - 1.1 x lO^Rg, - 0.15 
and ii - 105°, and at some stage they attain the values compatible with the S2-star parameters, i.e. a ^ 2.5 x 10^/?g and e ^ 0.87. 



Fine-tuning of the model parameters changes the quantitative 
results and allows us to decrease but these details cannot be 
compared with present-day observations. 

As explained above, jumps of mean eccentricity are caused 
by the Kozai process in an non-spherical (axially symmetric) 
perturbation of the central gravitational field. This we inter- 
pret as the influence of the disc, which simultaneously pro- 
vides the dissipative environment in a self-consistent manner. 
The process of Kozai-driven resonance prefers rather special 
orientation of initial trajectories, namely, almost perpendicular 
inclination. We can therefore speculate that the cuiTent fam- 
ily of S-stars was born jointly and these were originally mov- 
ing in a similar direction. They experienced the Kozai reso- 
nance and underwent the orbital decay ending in the present 
state of highly elliptic trajectories very near the black hole. 



Some of the nuclear stars in Sgr A* (Genzel et al. 120031 Levin 
& Beloborodov I2005t form a disc-like structure rather than a 
spherical distribution (nuclear stars are actually located in two 
rings, almost perpendicular to each other; Genzel et al. 2003 1. 
We can thus consider such rings as cuiTent prototypes of a for- 
merly existing ring in which S-stars were formed. 



4. Discussion 

In this paper we considered two complementary effects which 
both act as perturbations on otherwise free Keplerian motion, 
namely, gravity of the disc and the orbital decay due to cross- 
ing the disc slab. We computed a set of model trajectories of 
stars sinking towards a supermassive black hole in Sgr A*. We 
considered the region of gravitational dominance of the cen- 
tral black hole, however, we also took the gravity of the accre- 
tion disc and the corresponding dissipation into account. Stars 
do not strictly follow test-body (geodesic) motion in the cen- 
tral field, but instead they gradually sink to the centre and they 



experience occasional jumps of mean orbital elements. These 
second-order effects do not influence the motion on short time- 
scales of the order of one or several revolutions, nevertheless 
they are essential for long time-scales. Although the existence 
of a fossil disc in the Galaxy centre is a matter of debate, the 
idea offers a natural explanation for various pieces of observa- 
tional information, such as the occurrence of X-ray flares or the 
origin of young OB-type stars on highly elongated orbits. We 
therefore adopted the S2 orbital parameters in our simulations. 

A more complete picture will require not only a detailed 
knowledge of the disc structure but also coupling with other ef- 
fects, such as dynamical friction, gravitational radiation losses 
and the internal dissipation by tidal forces, which also con- 
tribute to the orbital decay, as well as the Lense-Thirring 
precession caused by rotation of the black hole. These and 
other mechanisms are beyond the scope of the present paper. 
Likewise, we neglected the effect of grazing encounters with 
other stars in the nuclear cluster All these processes will have 
to be taken into account in order to disentangle their impact 
on the evolving orbit and to measure black hole parameters in 
the future. Our model predicts that stars spend a considerable 
fraction of total time at the stage of large eccentricity. This can 
be reconciled with the current evidence which strongly prefers 
high eccentricities in the central arcsecond region of Sgr A*, 
because highly elliptical orbits are preferentially brought to 
small radii. 

The system of stars on energetically bound trajectories en- 
circling a supermassive black hole is a unique configuration 
allowing one to probe the gravitational field of the black hole 
by tracking the stellar proper motion. The available accuracy 
has justified the approximation of free motion, assuming that 
the stars behave as test particles that are insensitive to non- 
gravitational forces. This is fully satisfactory on the time-scale 
of several revolutions, i.e., of the order of tens to hundreds of 
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years. However, we showed that the effect of gaseous and/or 
dusty environment should not be neglected in the discussion 
of the long-term (Myr) evolution and that it can in fact help 
to reconcile the problem of the origin of the stellar population 
in the Galaxy centre. Although our approach involves various 
simplifications, the main message here is that one should take 
both effects into account simultaneously, i.e. the orbital dissipa- 
tion in the gaseous environment of the disc and its gravitational 
field. In fact, the idea of combining gravitational (conservative) 
and hydrodynamical (dissipative) star-disc interactions can be 
considered as a model case for a more general kind of process 
introducing a non-spherical perturbation to the central gravi- 
tational field and, simultaneously, a dissipative mechanism for 
the orbital decay. 
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